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The diffractive heavy quark cross sections are estimated considering photon-Pomeron (7-ff 5 ) inter- 
actions in hadron - hadron at RHIC, Tevatron, and CERN LHC energies. We assume the validity 
of the hard diffractive factorization and calculate the charm and bottom total cross sections and 
rapidity distributions using the diffractive parton distribution functions of the Pomeron obtained by 
the HI Collaboration at DESY-HERA. Moreover, we compare our predictions with those obtained 
using the dipole model. We verify that this process is a good test of the different mechanisms for 
diffractive heavy quark production at hadron - hadron colliders. 
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, A long-standing puzzle in the particle physics is the nature of the Pomeron (IP). This object, with the quantum 
ly-j 1 numbers of the vacuum, was introduced phenomenologically in the Regge theory as a simple moving pole in the 
complex angular momentum plane, to describe the high-energy behavior of the total and elastic cross-sections of the 
hadronic reactions Within the framework of the perturbative Quantum Chromodynamics (pQCD), the Pomeron 
is associated with the resummation of leading logarithms in s (center of mass energy squared) and at lowest order 
is described by the two-gluon exchange @. Due to its zero color charge the Pomeron is associated with diffractive 
events, characterized by the presence of large rapidity gaps in the hadronic final state. The diffractive processes have 
attracted much attention as a way of amplifying the physics programme at hadronic colliders, including searching for 
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New Physics (For recent reviews see, e.g. Refs. [1, The investigation of these reactions at high energies gives 
important information about the structure of hadrons and their interaction mechanisms. In particular, hard diffractive 
processes allow the study of the interplay of small- and large-distance dynamics within Quantum Chromodynamics 
(QCD). 

The diffractive processes can be classified as inclusive or exclusive events (See e.g. Q). In exclusive events, empty 
regions in pseudo-rapidity, called rapidity gaps, separate the intact very forward hadron from the central massive 
, object. Exclusivity means that nothing else is produced except the leading hadrons and the central object. The 
fT) ■ inclusive processes also exhibit rapidity gaps. However, they contain soft particles accompanying the production of a 
hard diffractive object, with the rapidity gaps becoming, in general, smaller than in the exclusive case. It is important 
to emphasize that at the LHC, the only possibility to detect double diffractive events is by tagging the intact hadrons 
in the final state, because of the pile up of events in each bunching crossing. It implies the key element to measure 
diffractive events at the LHC will be tagging the forward scattered incoming hadrons. 

A good testing ground for diffractive physics is the heavy quark production. In the last years, several authors 
have discussed this process considering different approaches. In particular, in Refs. [5|-|7| the total cross sections for 
Pomeron - hadron (inclusive single diffraction) and Pomeron - Pomeron (inclusive central diffraction) heavy quark 
production were calculated assuming the validity of the diffractive factorization formalism and that Pomeron has a 
partonic structure. The basic idea is that the hard scattering resolves the quark and gluon content in the Pomeron 8] 
and it can be obtained analysing the experimental data from diffractive deep inelastic scattering (DDIS) at HERA, 
providing us with the diffractive distributions of singlet quarks and gluons in the Pomeron Q . This model is usually 
denoted resolved Pomeron model. On the other hand, in Ref. [lfj the diffractive heavy quark production was 
estimated within the light-cone dipole approach, which is able to describe inclusive and diffractive processes in a 
unified way, with the key ingredient being the dipole-target cross section determined from HERA F% data. One of 
the main conclusions of that paper is the breakdown of factorization in diffractive hadronic collisions, which put in 
check the calculations presented in Refs. 0-0] • However, as both models describe the experimental data for single 
diffractive charm and bottom production, the correct approach to describe the inclusive diffraction process remains 
an open question. Finally, more recently, in Refs. (Til [l~2j . the exclusive double diffractive production of heavy 
quarks was estimated considering an improved version of the Durham model [rJ 03 (For a review see, e.g. fTo| ) and 
obtained that the cross section for this mechanism is much higher than the cross section for inclusive double-Pomeron 
interactions calculated in Consequently, the present scenario for diffractive heavy quark production is unclear, 
motivating the study of alternative processes which allow to constrain the correct description of the Pomeron. 

An alternative is the study of diffractive heavy quark production in hadronic interactions mediated by photons, 
usually denoted coherent or ultraperipheral collisions in the literature (For reviews see Ref. QJ]). Basically, in hadron 
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- hadron colliders, the relativistic hadrons give rise to strong electromagnetic fields. The photon stemming from the 
electromagnetic field of one of the two colliding hadrons can interact with one photon of the other hadron (two - 
photon process) or can interact directly with the other hadron (photon - hadron process). In the latter case, it is 
possible to study photon - parton (inclusive) and photon - Pomeron (diffractive) processes. Consequently, coherent 
processes allow us to study the inclusive and diffractive heavy quark photoproduction. In principle, the experimental 
signature of these processes is distinct and it can easily be separated. While in the inclusive heavy quark photon- 
hadron production we expect only one rapidity gap (associated to the photon exchange) and the dissociation of one 
of the incident hadrons, in diffractive heavy quark photoproduction, we expect the presence of two rapidity gaps 
and hadrons intact in the final state, similarly to exclusive double diffractive production. An important aspect to be 
emphasized is that in contrast to heavy quark production in hadron - Pomeron and Pomeron - Pomeron interactions, 
which have its predictions strongly affected by the soft absorption corrections, in photon - Pomeron interactions 
these corrections are expected to be very small. In Refs. [17H20l | heavy quark photoproduction was studied in detail, 
considering the dipole approach, and it was advocated that the cross sections for inclusive and diffractive interactions 
may be substantial (See also Refs. [HI, HH). I n these diffractive calculations the two-gluon perturbative model for 
the Pomeron was assumed. As the correct description of the Pomeron still is an open question, in this paper we 
will revisit this process and calculate the diffractive heavy photoproduction cross section assuming that the Pomeron 
has a partonic structure (resolved Pomeron model) and compare the predictions of both approaches for the total 
cross sections and rapidity distributions for hadronic collisions at RHIC, Tevatron and LHC energies. Our goal is 
twofold: (a) verify if the diffractive heavy quark photoproduction can be useful to discriminate between the resolved 
and perturbative models for the Pomeron and (b) compare our predictions with those obtained in Refs. pH Hoi ll^ . 

The paper is organized as follows. In the next section we present a brief review of coherent interactions in hadronic 
colliders. In Section [IIII we discuss the resolved Pomeron model and in Section [IVI we present our results. Finally, in 
Section [V] we summarize our main conclusions. 



II. COHERENT INTERACTIONS IN HADRONIC COLLIDERS 



Lets consider the hadron-hadron interaction at large impact parameter (b > Rf ll + Rh 2 ) an d at ultra relativistic 
energies. In this regime we expect the electromagnetic interaction to be dominant. In heavy ion colliders, the heavy 
nuclei give rise to strong electromagnetic fields due to the coherent action of all protons in the nucleus, which can 
interact with each other. In a similar way, it also occurs when considering ultra relativistic protons in pp(p) colliders. 
The photon stemming from the electromagnetic field of one of the two colliding hadrons can interact with one photon 
of the other hadron (two-photon process) or can interact directly with the other hadron (photon-hadron process). One 
has that the total cross section for a given process can be factorized in terms of the equivalent flux of photons of the 
hadron projectile and the photon-photon or photon-target production cross section [1 61 ] - In general, the cross sections 
for jh interactions are at least two order of magnitude larger than for 77 interactions (See e.g. Ref. 17]). In what 
follows our main focus shall be in photon - hadron processes. The cross section for the diffractive photoproduction of 
a final state QQ (= cc or bb) in a coherent hadron-hadron collision is given by, 
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is the equivalent photon flux, WA = zwy^NN 
hadron-hadron system. Considering the requirement that photoproduction is not accompanied by hadronic interaction 
(ultra-peripheral collision) an analytic approximation for the equivalent photon flux of a nuclei can be calculated, which 
is given by [l6j 
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where 77 = ujb/jL] Kq(ij) and Ki(rj) are the modified Bessel functions. Moreover, fj — uj (Rhi + Rh 2 )/lL and 
U(fj) — K\ (fj) — Kq (fj). The Eq. ([2]) will be used in our calculations of heavy quark production in pA collisions. On 
the other hand, for proton-proton collisions, we assume that the photon spectrum of a relativistic proton is given by 



dN 1 / p (u 

did 



2-kuj 



2uj 



VS. 



NN 



lnfi 



If 



1 



n 2n 2 3ft 3 



(3) 



3 



with the notation fl = 1 + [ (0.71 GeV 2 )/Q 2 lin ] and = ^/[^ (i _ 2 uj/JS^) } a (w/ 7 l) 2 . 

Some comments are in order here. Firstly, the coherence condition limits the photon virtuality to very low values, 
which implies that for most purposes, they can be considered as real. Moreover, if we consider pp/pPb collisions 
at LHC, the Lorentz factor is 7^ = 7455/4690, giving the maximum c.m.s. energy W 1P ~ 8390/1500 GeV. 
Therefore, while studies of photoproduction at HERA were limited to photon-proton center of mass energies of about 
200 GeV, photon-hadron interactions at LHC can reach one order of magnitude higher on energy. Secondly, in pA 
interactions one have that due to the asymmetry in the collision, with the ion being likely the photon emitter, the 
photon direction is known, which will implicate an asymmetry in the rapidity distribution (see below). Finally, in the 
resolved Pomeron model, in which the Pomeron has a partonic structure, the state X in the final state is associated 
to the hadronization of the Pomeron remnants. In contrast, in the perturbative Pomeron model analysed in [l9l l20l|. 
the final state X is not present. 



III. HEAVY QUARK PHOTOPRODUCTION IN THE RESOLVED POMERON MODEL 

In the resolved Pomeron model the diffractive cross sections are given in terms of parton distributions in the 
Pomeron and a Regge parametrization of the flux factor describing the Pomeron emission. The parton distributions 
has evolution given by the DGLAP evolution equations and arc determined from events with a rapidity gap or a 
intact proton, mainly at HERA. The basic idea to describe the diffractive heavy quark photoproduction in the resolved 
Pomeron model is that the cross section can be described similarly to the inclusive heavy quark photoproduction, with 
a diffractive gluon distribution g D (x,/j 2 ) replacing the standard inclusive gluon distribution (See e.g. Refs. [2~3, [25|). 
Explicitly one have 



(7(7/1 -> QQXh){W lh ) = / dxa^ w (W 7g )g n (x,^) (4) 

with Wjh the center-of-mass energy of the photon-hadron system, x the momentum fraction carried by the gluon, 
= Amq/W^ h and ttiq is the mass of the heavy quark. Moreover, the cross section for the photon-gluon fusion 



is given at leading order by [26j-|2 
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where eQ is the electric charge of the heavy quark, a em is the electromagnetic coupling constant and (3 — Armq^/W^ g . 
The distribution g D (x, /j, 2 ) can be obtained by the convolution of the Pomeron flux fp(xp) and the gluon distribution 
in the Pomeron gp(j3, /1 2 ): 

//* dx x 
dx P df3S(x - xpP)g P (P,n 2 )fp(xp) = / — —fp(xp)gp( , fj 2 ) , (6) 
J x Xp Xp 

where /i 2 = 4mg, xp is a fraction of the hadrons momentum transferred into the diffractive system and /3 is a fraction 
of the Pomeron momentum carried by the gluon participating in the diffractive scattering. The Pomeron flux fp(xp) 
enters in the form integrated over the four momentum transfer 

fp(xp) = J^ dtf p (x P ,t) , (7) 

where t m i n and t max are kinematic boundaries. In our calculations we will consider the Pomeron flux factor and gluon 
distribution in the Pomeron obtained by the HI Collaboration in its analyses of the diffractive structure function at 
HERA 9] , which assume a flux factor motivated by Regge theory with x p dependence given by, 

e B F t 

fp /p (xp,t) = A F ■ 2 , t) _ ± ■ (8) 



x 



The Pomeron trajectory is assumed to be linear, ap(t) = ajp(0) + a'pt, and the parameters Bp and a'p and 
their uncertainties are obtained from fits to HI FPS data 29]. Moreover, the Pomeron structure function has been 
modelled in terms of a light flavour singlet distribution £(z), consisting of u, d and s quarks and anti-quarks with 
u = d = s = u= d = s, and a gluon distribution g(z), with z being the longitudinal momentum fraction of the 
parton entering the hard sub-process with respect to the diffractive exchange, such that z — j3 for the lowest order 
quark-parton model process, whereas < (3 < z for higher order processes. 
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FIG. 1: Energy dependence of the rapidity distribution for the diffractive charm photoproduction in pp(p) collisions considering 
the resolved Pomeron model. 





RHIC 


TEVATRON 


LHC (7 TeV) 


LHC (14 TeV) 


u(pp — > p ® ccX <g> p 


47 nb 


214 nb 


709 nb 


1208 nb 



TABLE I: The total cross section for the diffractive charm photoproduction at the RHIC, TEVATRON and LHC energies 
considering the resolved Pomeron model. 



IV. RESULTS 



In what follows, we will compute the rapidity distribution and total cross sections for the diffractive photoproduction 
of heavy quarks in proton-proton and proton-nucleus at high energies. The resolved Pomeron model shortly reviewed 
in the previous section serve as input for the numerical calculations using Eq. (fTJ) for the energies of the current 
and future pp and pA accelerators. Namely, one considers \/Snn = 0.2(1.96) TeV for pp(p) collisions at RHIC 
(TEVATRON). At LHC we assume = 7 and 14 TeV for pp collisions and VSnn = 8.8 TeV for pPb collisions. 

Moreover, we assume in our calculations m c = 1.5 GeV and m& — 4.5 GeV. 

The distribution on rapidity Y of the produced open heavy quark state can be directly computed from Eq. |T]) , by 
using its relation with the photon energy cj, i.e. Y oc In (w/toq). Explicitly, the rapidity distribution is written down 
as, 

da [h 1 +h 2 ^h 1 ®QQX®h 2 )} dN l/hl {u) dN y/ha (u) 

dY = " du: <WQQfc 2 H + w doJ ^h^QQH, M ■ (9) 

where X is a hadronic final state resulting of the Pomeron fragmentation. 

In Fig. Q] we present our predictions for energy dependence of the rapidity distribution for diffractive charm 
photoproduction. As expected from the behaviour of the diffractive gluon distribution at small- a;, we have that the 
rapidity distribution increases with the energy. In Table U we present our predictions for the total cross section. 
In particular, we predict a value of 709 (1208) nb for pp collisions at 7 (14) TeV. Assuming the design luminosity 
C = 10' mb s -1 the corresponding event rates will be ~ 7000 (12000) events/second. In order to compare with 
previous estimates, in Fig. [5] we present our predictions for the diffractive charm and bottom photoproduction in 
pp collisions at y/s = 14 TeV (denoted resolved Pomeron) and those obtained in [20| using the dipole approach and 



the bCGC model [3l| for the dipole-proton cross section (denoted dipole-bCGC). For comparison we also present the 
prediction for the inclusive heavy quark photoproduction obtained using as input in the calculation the proton gluon 
distribution as given by the MRST parametrization [32]. Moreover, in Table [TTI we present our predictions for the 
total cross section. In the inclusive case, the values for charm and bottom production are either large, going from 
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FIG. 2: Rapidity distribution for the inclusive and diffractive charm (left panel) and bottom (right panel) photoproduction in 
pp collisions at \fs = 14 TeV. 







Inclusive 


Diffractive 
Dipole model 


Diffractive 
Resolved Pomeron 


Charm 


pp 


6697 nb 


161 nb 


1208 nb 




pPb 


5203 fib 


145 fib 


694 fib 


Bottom 


pp 


123 nb 


0.52 nb 


15 nb 




pPb 


55 fib 


0.2 fib 


4.5 fib 



TABLE II: Comparison between the total cross sections for the inclusive and diffractive charm and bottom photoproduction 
in pp and pPb collisions at LHC. 



some units of nanobarns at bottom to microbarns at charm. Therefore, these reactions can have high rates at the 
LHC kinematical regime. On the other hand, the cross sections for diffractive production considering the resolved 
Pomeron model are a factor 5 smaller than the inclusive case. In comparison with the predictions of the dipole 
model, our predictions are a factor > 7 larger. It is important to emphasize that although the smaller values of the 
cross sections for diffractive production in comparison to the inclusive case, due the clear experimental signature of 
the diffractive processes, two rapidity gaps and two intact hadrons in the final state, the experimental analyses is 
expected to be feasible. 

The diffractive heavy quark photoproduction can also be studied in proton - nucleus collisions at LHC, which are 
expected to occur in the next year. In this case the cross sections are enhanced by a factor Z 2 present in the photon 
flux of a nuclei. In Fig. [3] we present our predictions for the inclusive and diffractive heavy quark photoproduction 
considering pPb collisions at -y/i = 8.8 TeV. In this case the rapidity distribution is asymmetric due to the dominance 
of the ion as the photon emitter. The large difference between the predictions from different models observed in pp 
collisions also is present in the pA case. In particular, as observed in Table lU we predict large values for the total 
cross section, which becomes the experimental test of our results feasible. It is important to emphasize that the 
predictions of the dipole model for pA collisions are presented here for the first time. Finally, we would like to point 
out that the diffractive photoproduction should also occur in AA collisions, but in order to estimate the cross section 
using the resolved Pomeron model, we would need to specify the diffractive parton distribution of the Pomeron and 
its flux in a nuclei, which are currently not constrained by experimental data. We postpone the analysis of this case 
for a future publication. 

Lets discuss the experimental separation of the diffractive photoproduction of heavy quarks and compare our 
predictions with those obtained considering Pomeron - Pomeron interactions . In comparison to the inclusive heavy 
quark hadroproduction (See e.g. [HI), which is characterized by the process p + p —> X + QQ + Y, with both 
proton producing hadronic final states, the photoproduction cross sections are a factor m 10 3 smaller. However, as 
diffractive photoproduction is an exclusive reaction, p + p ^ p® QQ ® p, the separation of the signal from hadronic 
background would be relatively clear. Namely, the characteristic features in photoproduction at coherent collisions 
are low pr heavy quark spectra and a double rapidity gap pattern. Moreover, the detection (Roman pots) of the 
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FIG. 3: Rapidity distribution for the inclusive and diffractive charm (left panel) and bottom (right panel) photoproduction in 
pPb collisions at -Js = 8.8 TeV. 

scattered protons can be an additional useful feature. In hadroproduction, the spectra on transverse momentum of 
heavy quarks are often peaked around heavy quark pair mass, pt ~ m QQ- Therefore, we expect that a cut in the 
transverse momentum of the pair, for instance pr < 1 GeV, could eliminate most part of the contribution associated 
to the hadroproduction of heavy quarks (See e.g Refs. EH). Furthermore, the rapidity cut would enter as an 
auxiliary separation mechanism. This procedure is specially powerful, since there will be a rapidity gap on both 
sides of the produced heavy quarks in the diffractive case. Moreover, in comparison with the hadroproduction of 
heavy quarks, the event multiplicity for photoproduction interactions is lower, which implies that it may be used as 
a separation factor between these processes. 

In the case of diffractive photoproduction of heavy quarks we expect the presence of two rapidity gaps in the final 
state, similarly to two - photon or Pomeron - Pomeron interactions. Consequently, it is important to determine the 
magnitude of this cross section in order to estimate the background for these other channels. Lets start considering 
proton - proton collisions. In comparison to production of heavy quark in two - photon interactions [34| . our predictions 
are a factor m 10 4 larger, which implies that this contribution can be disregarded. In contrast, in comparison to the 
inclusive double diffractive heavy quark production studied in Ref. [f|, which predict <j P g p ~ 18 /ib and cr^ p ~ 0.5 
/jb for LHC energies, our predictions are a factor > 10 smaller. In comparison to the exclusive double diffractive 
production of charm quarks presented in [l2| at Tevatron energy, it is a factor 30 smaller. It is important to 
emphasize that these predictions are strongly dependent on the value used for the gap survival factor, while our 
results should not be modified by soft absorption corrections. 

In the case of proton - nucleus collisions, the magnitude of the heavy quark production considering Pomeron - 
Pomeron interactions still is an open question. In particular, the value of the gap survival factor in nuclear collisions 
is a subject that deserves a detailed study. Preliminary studies, as in Ref. |35|, indicate that this factor is very small 
(«8x 10~ 4 ), which implies that jlP interactions can dominate the heavy quark production in pA collisions. 

Finally, although our predictions are smaller than those obtained considering double pomeron interactions, it is 
expected that emerging hadrons from double Pomeron processes have a much larger transverse momentum than those 
resulting from diffractive photoproduction processes. Consequently, in principle it is possible to introduce a selection 
criteria to separate these two processes. However, this subject deserves more detailed studies. 

V. CONCLUSIONS 

In summary, we have computed the cross sections for diffractive photoproduction of heavy quarks in pp and pA 
collisions at LHC energies. This has been performed using the resolved Pomeron model based on the diffractive factor- 
ization formalism and that the Pomeron has a partonic structure, which describe quite well the HERA experimental 
data for the diffractive structure function. The obtained values are shown to be sizeable at LHC. The feasibility of 
detection of these reactions is encouraging, since their experimental signature should be suitably clear. Furthermore, 
considering the large difference between the predictions of the distinct models, they enable to constrain the underlying 
model for the Pomeron, which is fundamental to predict the observables which will be measured in hadron-hadron 
collisions at LHC. 
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